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Abstract. We present new results for Casimir forces between rigid bodies which 
impose Dirichlet boundary conditions on a fluctuating scalar field. As a universal 
computational tool, we employ worldline numerics which builds on a combination of the 
string-inspired worldline approach with Monte-Carlo techniques. Worldline numerics 
is not only particularly powerful for inhomogeneous background configurations such 
as involved Casimir geometries, it also provides for an intuitive picture of quantum- 
fluctuation-induced phenomena. Results for the Casimir geometries of a sphere above 
a plate and a new perpendicular-plates configuration are presented. 



1. Introduction 

Casimir energies and forces are geometry dependent. Determining the geometry 
dependence is a challenge both experimentally as well as theoretically. Computing the 
geometry dependence of Casimir energies can be viewed as a special case of the more 
general problem of evaluating the effects of quantum fluctuations in a background field 
V{x). For instance, the space(-time) dependence of the background field can then be 
used to model the Casimir geometry. 

A universal tool to deal with quantum fluctuations in background fields is given by 
the effective action T[V] which is the generating functional for IPI correlation functions 
for V. In the present work, we consider a fluctuating real scalar quantum field (p 
interacting with the background potential according to ~ V{x)(j)'^. For this system, 
the effective action can be evaluated from 



Here h and c are set to 1. The integral over the Gauf^ian fluctuations boils down to 
a sum over the spectrum {A} of quantum fluctuations. This spectrum consists of the 
eigenvalues of the fluctuation operator. 



The relation to Casimir energies becomes most obvious by confining ourselves to time- 
independent potentials ^(x). Then, the sum over the time-like component of the 
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spectrum can be performed. In Euclidean spacetime, we use —d"^ = —d^—V^, —df pf, 
and the summation/integration over pt results in 

^M-^ = ^E-. (3) 

where = \^ —pi denote the spatial (pt-independent) part of the fluctuation spectrum. 
Here, we have defined the Casimir energy from the effective action by dividing out the 
extent Lt of the Euclidean time direction. The relation to a sum over "ground-state 
energies" ^ \hjj now becomes obvious. 

The general strategy for computing E\y] seems straightforward: determine the 
spectrum of quantum fluctuations and sum over the spectrum. However, this recipe 
is plagued by a number of profound problems: first, an analytic determination of the 
spectrum is possible only in very rare, mainly separable, cases. Second, a numerical 
determination of the spectrum is generally hopeless, since the spectrum can consist of 
discrete as well as continuous parts and is generically not bounded. Third, the sum over 
the spectrum is generally divergent, and thus regularization is required; particularly in 
numerical approaches, regularization can lead to severe stability problems. And finally, 
an unambiguous renormalization has to be performed, such that the physical parameters 
are uniquely fixed. 

For a solution of these problems, the technique of worldline numerics has been 
developed P and has first been applied to Casimir systems in j2]. This technique is 
based on the string- inspired approach to quantum field theory [Sj. In this formulation, 
the (quantum mechanical) problem of finding and summing the spectrum of an operator 
is mapped onto a Feynman path integral over closed worldlines. For the present scalar 
case, the effective action then reads 

r|V| = -i r f e-'-'-A/- 

Now, the effective action, or, more specifically, the Casimir energy, is obtained from 
an integral over an ensemble of closed worldlines in the given background V{x). This 
seemingly formal representation can be interpreted in an intuitive manner: A worldline 
can be viewed as the spacetime trajectory of a quantum fluctuation. The auxiliary 
integration parameter T is called the propertime and specifies a fictitious "time" which 
the fluctuating particle has at its disposal for traveling along the full trajectory. Larger 
values of the propertime thus correspond to worldlines with a larger extent in spacetime; 
hence, the propertime also corresponds to a smooth regulator scale, with, e.g., short 
propertimes being related to small- distance UV fluctuations.! 

Most importantly from a technical viewpoint, the problem of finding and summing 
over the spectrum is replaced by one single step, namely taking the path integral. 
Moreover, for any given value of propertime T this path integral is finite. Possible 
UV divergencies can be analyzed with purely analytical means by studying the small-T 

X For reasons of definiteness, wc have therefore cut off the propertime integral at the lower bound at 
l/A^ with the UV cutoff scale A. 



'-+V[x{t))) 
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behavior of the propertime integral, and thus no numerical instabilities are introduced 
by the regularization procedure. The renormalization can be performed in the standard 
way, for instance, by analyzing the corresponding Feynman diagrams with the same 
propertime regulator and by fixing the counterterms with the aid of renormalization 
conditions accordingly. 

Further advanced methods which can deal with involved Casimir configurations 
have been developed during the past years, each with its own respective merits. In 
particular, we would like to mention the semiclassical approximation |3], a functional- 
integral approach using boundary auxiliary fields and the optical approximation [H]. 
These methods are especially useful for analyzing particular geometries by purely or 
partly analytical means. 

Of course, a purely analytical evaluation of the path integral again is only possible 
for very rare cases, but a numerical evaluation is straightforwardly possible with Monte 
Carlo techniques and can be realized with conceptually simple algorithms, as described 
in the next section. Applications to Casimir geometries will be presented in Sect. Eland 
conclusions are given in Sec. |3] 

2. Worldline numerics for Casimir systems 

With the aid of the normalization of the path integral jj, we note that (0} can be 
written as 

where the subtraction of —1 ensures that r[y = 0] = 0. The expectation value in 
has to be taken with respect to the worldline ensemble, 

Vx ...e--J^''^A ( [ Vxe-yo'^-A''. (6) 

, -(T)=x(0) J \Jx{T)=x{0) J 

In the present work, we focus on the "ideal" Casimir effect induced by real scalar 
field fluctuations obeying Dirichlet boundary conditions; i.e., the boundary conditions 
are satisfied at infinitely thin surfaces. This situation can be modeled by choosing 
V{x) = g Jy; da S^^\x — Xa), where da denotes the integration measure over the surface 
S, with Xcr being a vector pointing onto the surface. The Dirichlet boundary condition 
is then strictly imposed by sending the coupling g to infinity, g ^ oo [HIH]. 

Moreover, we are finally aiming at Casimir forces between disconnected rigid 
surfaces, which can be derived from the Casimir interaction energy, 

^Casimir = E[Vi + V2] " E[V^] - £[¥2], (7) 

where we subtract the Casimir energies of the single surfaces Vi and V2 from that of 
the combined configuration Vi + V2;§ the former do not contribute to the force. This 

§ Alternatively, the single-surface subtraction terms in (O can be viewed as subtracting the Casimir 
energy for the surfaces at infinite separation. 
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definition, together with the Dirichlet boundary condition, leads to 



^Casimir = " 2 / ^ ^ {^v[A) x ' 




(8) 



where 0v[x] = 1 if a given worldhne intersects both surfaces E = Ei + E2 represented 
by the background potentials V = Vi + V2, and 6y[x] = otherwise. This recipe has 
a simple interpretation: any worldline which intersects both surfaces corresponds to 
a quantum fluctuation that violates the Dirichlet boundary conditions. Its "removal" 
from the set of all fluctuations contributes "one unit" to the negative Casimir interaction 
energy. 

The worldline numerical algorithm is based on an approximation of the path integral 
by a sum over a finite ensemble of tt-l number of paths, each of which is characterized 
by N discrete points per loop (ppl). These points are obtained by a discretization of 
the propertime parameter on each loop, Xj = x(rj), i = 1,...,N, with (xj)^ G R. 
For an efficient generation of the worldline ensemble which obeys a Gauf^ian velocity 
distribution required for various algorithms are available, see [SEl- 

In summary, worldline numerics offers a number of advantages: first, the whole 
algorithm is independent of the background; no particular symmetry is required. Second, 
the numerical cost scales only linearly with the parameters riL, A^, and the dimensionality 
of the problem. The numerically most expensive part of the calculation is a diagnostic 
routine that detects whether a given worldline intersects both surfaces or not, returning 
the value 0v[x] = 1 or 0, respectively. Optimizing this diagnostic routine for a 
given geometry can lead to a significant reduction of numerical costs. Details of this 
optimization for the geometries considered below will be given elsewhere [Til] . 

3. Application to Casimir geometries 
3.1. Sphere above plate 

The geometry of a sphere above a plate is the most relevant configuration as far as recent 
and current experiments are concerned Therefore, also worldline numerics has first 
been applied to this case |2]. Here we extend these studies, arriving at significantly 
improved results with much smaller error bars and for a wider range of parameters. In 
the following, we exclusively discuss the massless case, m = 0. 

It is interesting to compare our results to the proximity force approximation (PEA) 
[T^ which is the standard tool for estimating the effects of departure from planar 
geometry for Casimir effects. In this approach, the curved surfaces are viewed as a 
superposition of infinitesimal parallel plates, and the interaction energy is then obtained 



Here, EppA denotes a "suitable" auxiliary surface in between the Casimir surfaces Ei 
and E2, and da is the corresponding surface element of EppA- The distance d between 



by 
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two points on Si and S2 has to be measured along the normal to SppA- Obviously, the 
definition of -EppA is ambiguous, owing to possible different choices of SppA- The two 
extreme cases are SppA = Si or SppA = S2. The difference in £'pfa for these two cases 
is considered to represent a rough error estimate of the PFA. In the above formula, E-p-p 
denotes the classic parallel-plate result for the energy per unit area A 

with a denoting the plate separations, and cpp = 2 for an electromagnetic (EM) field or 
a complex scalar, and cpp = 1 for the present case of a real scalar field fiuctuation. 

For the configuration of a sphere above a plate, we can choose SppA equal to the 
plate (plate-based PFA), or equal to the sphere (sphere-based PFA) as the extreme 
cases. The corresponding results of the integral of can be given in closed form, see, 
e.g., [H]. In the limit of small distances a compared to the sphere radius i?, a/R <^ 1, 
both PFA's agree, 

EpFA(a/i?«l) = -cpp^^. (11) 

It is useful to display the resulting Casimir energies normalized to this zeroth-order 
small-separation PFA limit as it is done in figure ^ The dashed and dot-dashed lines 
depict the plate-based and sphere-based PFA, respectively. For larger separations, both 
cases predict a decrease of the Casimir energy relative to the zeroth-order PFA in (|TT|l . 

Our numerical worldline estimate confirms the zeroth-order PFA in the limit of 
small a/R. But in contrast to the PFA curves, worldline numerics predicts a relative 
increase of the Casimir energy in comparison with for larger separations/smaller 
spheres. We conclude that the PFA should not at all be trusted beyond the zeroth 
order: the first-order correction does not even have the correct sign. As a most 
conservative estimate, we observe that the PFA deviates from our result by at least 
1% for a/R > 0.01, which confirms and strengthens the result of A more detailed 
investigation of the validity bounds of the PFA will be given elsewhere ^Oj- It is 
interesting to observe that the zeroth-order PFA (fTT|) still seems to be a reasonable 
estimate up to a/R ~ 0.1, indicating that the true curvature effects compensate for the 
higher-order PFA corrections. 

Finally, we note that our results agree with the optical approximation jB] for 
a/R < 0.1, confirming the absence of diffractive effects in this regime, which are 
neglected by the optical approximation. For even larger separations a, we observe a 
monotonous increase of the Casimir energy relative to (jllj) . In this regime, our results 
agree quantitatively with those obtained from the "KKR" multi-scattering map method 
presented by A. Wirzba at this QFEXT05 workshop ^1]. Most importantly, we do not 
observe a Casimir-Polder law for large a/R, which would manifest itself in an {a/R)~'^ 
decrease in figure ^ at the large- (a/i?) side. Since a Casimir-Polder law is expected for 
the electromagnetic case, our results for the Dirichlet scalar provide clear evidence for 
the fact that the relation between Casimir forces for the EM field and for the Dirichlet 
scalar is strongly geometry dependent. In fact, the latest results of Ref. [14] include an 
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Figure 1. Casimir energy for the sphere-plate configuration normalized to the zcroth- 
order PFA formula Hll|l : the dashed and dot-dashed lines depict the plate-based and 
the sphere-based PFA estimates, respectively. The circle symbols display our worldline 
numerical result. The deviation from the PFA estimate characterize the relevance of 
Casimir curvature effects. Also shown is the result from the optical approximation 
[S], which, within its validity limits a/ R < 0.1, agrees well with our result. For larger 
a/i?, we find satisfactory agreement with the "KKR" multi-scattering map method 
presented at this workshop [Tl). 



analytic proof that the energy ratio of Fig. ^approaches a constant, ISO/tt^ ~ 1.848 
for large a/ R for the Dirichlet scalar. The analysis of the sphere-plate configuration for 
the EM field therefore still remains an open unsolved problem. Note that this does not 
affect our conclusions about the PFA, since also the PFA does not treat the EM case or 
the Dirichlet scalar in a different manner. 

3.2. Perpendicular plates 

A particularly inspiring geometry is given by a variant of the classic parallel-plate case: 
a semi-infinite plate perpendicular to an infinite plate, such that the edge of the semi- 
infinite plate has a minimal distance a to the infinite plate, see figure El (left panel). 
Whereas Casimir's parallel-plate case has only one nontrivial direction (the one normal 
to the plates), this perpendicular-plates case has two nontrivial directions but still only 
one dimensionful scale a. This fixes the scale dependence of the energy per unit length 
unambiguously, 

-T^ = -^^1440^' ^^^^ 
where Lt denotes the extent of the system along the remaining trivial transversal 
direction. The unknown coefficient 7^ results from the effect of quantum fluctuations 
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Figure 2. Left panel: sketch of the perpendicular-plates configuration with (an artist's 
view of) a typical worldline that intersects both plates. Right panel: density plot 
of the effective action density C for the perpendicular plates case; the plot shows 
ln(2(47ra^)^|>C|). The position of the perpendicular plates are indicated by solid lines 
for illustration. 



in this geometry and will be determined by worldline numerics. 

Let us first note that the PFA does not appear to be useful for the perpendicular- 
plates case, because the surfaces cannot reasonably be subdivided into infinitesimal 
surface elements facing each other from plate to plate. For instance, choosing either of 
the plates as the integration surface in 0, the PFA would give a zero result. However, 
in the worldline picture, it is immediately clear that the interaction energy is nonzero, 
because there are many worldlines which intersect both plates, as sketched in figure 
El (left panel). As a direct evidence, we plot the negative effective action density C 
(effective Lagrangian) in figure El (right panel); the effective action is obtained from 
r = J di^xC Brighter areas denote a higher density of the center-of- masses of those 
worldlines which intersect both plates. 

Integrating over the effective action density, we obtain the universal coefficient 

7^ = 0.87511 ±0.00326, (13) 
using ul = 40000 worldlines with = 200000 ppl generated by the v loop algorithm 

m 

4. Conclusions 

We have presented new results for interaction Casimir energies, giving rise to Casimir 
forces between rigid bodies, induced by a fluctuating real scalar field that obeys Dirichlet 
boundary conditions. We have used worldline numerics as a universal tool for dealing 
with quantum fluctuations in inhomogeneous backgrounds. 
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For the experimentally relevant sphere-plate configuration, we have performed 
extensive numerical studies, confirming earlier findings [2] with a significantly higher 
precision and narrowing the validity bounds of the proximity force approximation even 
further. Moreover, our results for small spheres for the Dirichlet scalar shows no sign 
of a Casimir-Polder law, as it would be expected for the EM field. This provides clear 
evidence for a different role of Casimir curvature effects for these two different field 
theories, leaving the sphere-plate configuration with a fiuctuating EM field as a pressing 
open problem. 

Furthermore, we have investigated a new geometry of two perpendicular plates 
which has been inaccessible so far for other approximation techniques. The configuration 
is representative for a whole new class of Casimir systems involving sharp edges, where 
diffractive portions of the fiuctuating field will play a major role. 
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